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We study the bound states of two spin-| fermions interacting via a contact attraction (character- 
ized by a scattering length) in the singlet channel in 3D space in presence of a uniform non-Abelian 
gauge field. The configuration of the gauge field that generates a Rashba type spin-orbit interaction 
is described by three coupling parameters {\x,\y,\z). For a generic gauge field configuration, the 
critical scattering length required for the formation of a bound state is negative, i.e., shifts to the 
"BCS side" of the resonance. Interestingly, we find that there are special high-symmetry configura- 
tions (e.g., = Ay = Az) for which there is a two body bound state for any scattering length however 
small and negative. Remarkably, the bound state wave functions obtained for high-symmetry con- 
figurations have nematic spin structure similar to those found in liquid '^He. Our results show that 
the BCS-BEC crossover is drastically affected by the presence of a non-Abelian gauge field. We 
discuss possible experimental signatures of our findings both at high and low temperatures. 
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I. INTRODUCTION 

Quantum emulation experiments with cold quantum 
gasei^l^ and optical lattices hold the promise of provid- 
ing clues to understanding many outstanding issues of 
quantum condensed matter physics such as high temper- 
ature superconductivity, quantum hall effect, etc, and the 
high energy physics of strongly coupled gauge theories.!^ 
While this has led to a flurry of activity, many exper- 
imental challenges remain in the way of redemption of 
this promise. Particular among them are the problem 
of entropy removal and the creation of magnetic (gauge) 
fields. 

Realization of magnetic fields has been achieved by 
rotation^, however, attaining magnetic fields correspond- 
ing to quantum hall regimes has serious experimen- 
tal challenges. There have been many theoretical sug- 
gestions for the generation of artificial gauge fielda^J'^, 
both abe lian a nd non-abelian. Recently Spielman and 
coworker^iSIlil used Raman coupling between hyperfine 
states to produce synthetic gauge fields. They studied 
Bose condensates of ^^Rb atoms and investigated the 
punching in of vortices when a U(l) gauge field cor- 
responding to a magnetic field is tuned. Depending 
on the degeneracy of the lowest Raman coupled states, 
one can also generate non-abelian gauge fields. The 
condensatio n of b osons in non-abelian fields have been 
investigatedi ^^ l ^^ l 

These developments provide us the motivation to study 
fermions in non-abelian gauge fields. The simplest in this 
class is the case of spin-^ particles coupled to an SU(2) 
gauge field. Study of such systems within the cold atoms 
context will enable experimental realization and under- 
standing of fermionic Hamiltonians with spin orbit inter- 
actions that can lead to interesting topological phases of 
matter.lSMl 

Readers who wish to obtain a qualitative understand- 
ing of our work may read Sec. ITT] where we state our 



problem and summarize our results, followed by Sec. [V 
which discusses the significance of these results. Sec. [HI 



contains details of our calculations, and Sec. |IV| provides 
a qualitative discussion of the physics of our results. 



II. STATEMENT OF THE PROBLEM AND 
SUMMARY OF RESULTS 

We consider spin-^ fermions moving in 3D continuum 
in a non-abelian gauge field. The simplest realization of 
this is described by the Hamiltonian 
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GF 



*(r) 
(1) 



where ^'(r) — {ip„{r)},a =tii is a two component 
spinor field (spin quantization along z-axis), pi is the 
momentum operator (i = x,y, z), 1 is the SU(2) identity, 
are Pauli spin operators (/i = x,y,z), describe 
a uniform gauge field. We work with units where the 
mass of the fermions and h are unity. Indeed even a uni- 
form non-Abelian field leads to interesting physicJ^^l^ an 
example of which we demonstrate in this paper. 

Motivated by the recent experiments mentioned above, 
we consider non-Abelian gauge fields of the type A'^ — 
Ai(5f which leads to a generalized Rashba Hamiltonian 
describing an anisotropic spin orbit interaction 



■Hr = Jd'r^ir) *(r) 

PX = KPxSx + KjPySy + XzPz^z- 



(2) 



Here an inconsequential constant term has been dropped. 
The gauge coupling strength is A = y^A^ + A^ + A^, and 

the vector A = A A = A^jfi-c + \yey + Xz^z defines a gauge 
field configuration (GFC). 
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FIG. 1. (color online) (a) Two-body phase diagram in the gauge-field configuration (GFC) space described in (|2|. Lines in 
dark blue (along the axes) correspond to extreme prolate (EP) GFCs, dark green (along 45° to the axes on the coordinate 
planes) correspond to extreme oblate (EO) GFCs, and the red line (along the body diagonal) corresponds to spherical (S) 
GFC. The wings in violet hue correspond to oblate GFCs. For the oblate, including EO GFCs, and the spherical (S) GFC a 
bound state is obtained for any scattering length i.e., a^c = 0~. Regions with a blue hue, including the axes (EP), correspond 
to GFCs that require a non-zero critical scattering length for the formation of a bound state, (b) The critical scattering length 
along the path EP-S-EO as shown in (a). For the EP GFCs, the symmetry of the bound state wave function corresponds 
to an extended Balian-Werthamer (BW) state with a biaxial nematic spin order, while that in the EO state corresponds to 
an extended Anderson-Brinkman-Morel (ABM) state with a uniaxial nematic spin order. The state evolves smoothly from a 
biaxial nematic to a uniaxial nematic passing through the spherical configuration (S) where the bound state is rotationally 
symmetric. 



We now describe the interaction between fermions by 
a contact attraction model in the singlet channeP^ 

= ^1 dV5t(r)5(r) (3) 

where S^{r) is the singlet creation operator, and v is 
the bare contact interaction. The theory described by 
the Hamiltonian "H = Hr + Hv requires an ultraviolet 
momentum cut-off A. The bare contact interaction pa- 
rameter V is A-dependent and satisfies the regularization 
relation - 4- A = . ^ , where a, is the s-wave scattering 
length in the absence of the gauge field ("free vacuum"). 
It is well knowrfi^ that for a pair of spin- ^ fermions in 
free vacuum, there is no bound state when < (con- 
ventionally called the "BCS side"), and a bound state 
develops when I /as — ("resonance"), and for > 
("BEC side"), the binding energy Eh — This result 

embodies the fact that a critical attraction, characterized 
by the critical scattering length Asc, is needed to obtain 
a two-body bound state in the 3D free vacuum where 
l/asc = 0. 

In this paper we address the question of how a uni- 
form gauge field described by ^ affects the nature of the 
bound state of two fermions interacting via ^. To this 
end, we obtain the "phase-diagram" of the two-fermion 
problem in the GFC space described by the parameters 
Xx,\y,Xz of ([2]), by studying the bound state as a func- 
tion of the free vacuum scattering length . 

GFCs can be conveniently classified as being prolate 



when two of the As are equal and smaller than the third, 
oblate when two of the As are equal and larger than the 
third, spherical (S), when all three As are equal, and 
generic when no two As are equal. Our main findings 
are summarized in Fig. [T] We show that for prolate and 
generic GFCs, the critical scattering length a^c required 
for the bound state formation is negative i.e., shifts to 
the BCS side. However, for oblate and spherical GFCs 
Use vanishes, i.e., there is a bound state for any scatter- 
ing length (see Fig.[ljb)). The key difference between the 
oblate and spherical cases is the size of the binding en- 
ergy of the bound state. In the deep BCS side, for oblate 
gauge fields, the binding energy depends exponentially 
on fls and A, while for spherical gauge fields, an alge- 
braic dependence is obtained. Evidently, these results of 
the two-body problem suggest that many body physics 
of fermions, in particular, the crossover from the BCS 
regime to the BEC regime will be spectacularly affected 
by the presence of a non-abelian gauge field. Moreover, 
our results below indicate that the superfiuid obtained at 
low temperatures will also have additional spin nematic 
order induced by the gauge field. 



III. TWO-BODY PROBLEM IN PRESENCE OF 
NON-ABELIAN GAUGE FIELDS 

For any GFC, the single particle states are described 
by the quantum numbers of momentum k and helicity a 
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(which takes on values ±): 



\ka) = |fe) (8) \akx) 



(4) 



where |fc) is the usual plane- wave state, and \ak\) is the 
spin coherent state in the direction ak\ , with k\ defined 
analogous to pa of ([2|. The two helicity bands disperse 



fc2 
2 



a\k\ 



(5) 



The full two body Hamiltonian % generically has only 
two symmetries - global translation and time reversal. 
Therefore, the only good quantum number is the center 
of mass momentum of the two particles. We shall fo- 
cus attention on states with zero center of mass momen- 
tum, and perform a T-matrix analysiJ^^ in the relative- 
momentum and helicity bases. The components of the 
T matrix have the matrix structure Tg^'(w), where lo is 

the energy, and run over (++,H — , — h, ), the 

helicity indices of the two fermions. Since the interaction 
is only in the singlet channel, it follows that components 
with indices (H — , — h) vanish. This analysis, along with 
the regularization discussed earlier, readily provides the 
condition for bound-state formation 
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where V is the volume of space under consideration. Iso- 
lated poles of the T matrix, which correspond to bound 
states, are obtained by finding the roots £^ of (|6]). We 
shall now present results for particular GFCs including 
the nature of the bound-state wave functions. 



A. Extreme prolate (EP) 

In EP GFCs, two of the gauge couplings vanish (say 
Ax = Ay = 0) while only one is nonzero (Az = A). 
Such configurations correspond to the axes marked in 
blue (along the axes) in Fig. [IJa). These GFCs pos- 
sess, in addition to translation and time reversal, spatial 
and spin rotation symmetries about the z axis. The one 
particle dispersion ([s]), for this case, provides the scatter- 
ing threshold Eth = — A^. Defining the binding energy 
Eb = —{E — Eth), we find from the solution of (|6| that a 
bound state appears only for positive scattering lengths 
(Fig. [2]), with 



Eh = 



> 0. 



(7) 



The critical scattering length corresponds to resonance 
i.e., l/flsc = 0. 

These results for Ei, and are identical to those of the 
two-body problem in free vacuum. There is, however, a 
crucial difference. The wave function of the bound state 
in the absence of the gauge field is a spin singlet. In 



an extreme prolate gauge field, the bound state wave 
function has two pieces, 

^b^Mr)\n-it) + (r ) in + ;t> , (s) 

where Mr\ = Eka and Mr) - Efc„ 

are, respectively, symmetric and anti-symmetric func- 
tions of the relative coordinate r. The first piece is 
a spin singlet, while the second piece (which vanishes 
when A — >■ 0), has a triplet spin wave function. This 
wave function corresponds to an extended BW stat^^ 
of the B-phase of ^He with an additional singlet piece. 
This state has a spin-nematic ordei^ corresponding to 
a biaxial nematic, consistent with the symmetries of the 
Hamiltonian. 



B. Extreme oblate (EO) 

These configurations have one of the gauge couplings 
equal to zero, and the other two equal and non-zero, and 
are marked by green lines (lines at 45° to the axes on 
the coordinate planes) in Fig. [T] We consider the case 
with Ax = Aj, = Az = 0. For these GFCs, we have in 
addition to translation and time reversal, a symmetry of 
global (spatial -|- spin) rotation about the z axis gener- 
ated hy Jz = Lz + -^Tz, where Lz is the z component of 
the orbital angular momentum operator. 

The secular equation Q, in this case, reduces to 



V2 
Aa., 



2Eb 
A2 




' - log I — — 
2 '^\2Eb 



The results presented for the EO case in the published 
version of this paper (Phys. Rev. B 83, 094515 (2011)) is 
based on the analysis of the above equation. The analy- 
sis, however, was performed with an erroneous factor of 
2 in the last term: —^^oe{^) was used instead of the 
correct term — ^log(^|^). The results shown below (in- 
cluding Fig. [2| are the correct results for the EO GFCs. 
We note that no qualitative conclusions and physics dis- 
cussed in the published version is altered by this error, 
and in addition, the results presented in the published 
version for all other GFCs are correct. 

The solution of ^ provides an interesting result: there 
is a bound state for any scattering length, negative or 
positive, i.e., a^c = 0~. For small negative (BCS 
regime), we obtain the binding energy (referred to the 
scattering threshold Etji — —4-) ^ 



Eb 
A2 



n=0 



-2^/2 



2 -2n/2 



(9) 



where a„s are positive rationals tending asymptotically 
to (4e)5; ag = l;ai = l;a2 = j;a3 — ^. Thus the 
binding energy is exponentially small for small negative 
flj. For small positive (BEG regime), we have 

Eb ^ 1 ^ 1 _ (Aa,)2 
A2 



12 



(10) 
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which recovers the binding energy of 1 /a^ in the hmit of 
free vacuum (A — )■ 0). When l/a^ = (resonance), the 
binding energy is determined solely by A; we obtain, near 
resonance. 



A2 



C 



2V2C 1 



4C(1+C) 1 
Vl + 2CXas ^ (1 + 2C)2 [XasY 



.(11) 



where C ~ 0.2196. The full evolution of the bound state 
energy as a function of the scattering length Os is shown 
in Fig.[2| 

The bound-state wave function, again, has two pieces 

*b (X v.(r)i n - ;t> + Mr)\ n) + ^:(r-)i u) (12) 

where, 7/;,(r) = - Efca 2eT„-£ ' ^^^^ Va(r) - 

with the a; axis. The first piece is orbitally symmetric 
{tpsif)) spin singlet (the first term), and the second piece 
(next two terms) consists of an antisymmetric orbital 
wave function (ipair)) and a spin structure correspond- 
ing to that of the ABM statd^l in the A-phase of ^He. 
This state has uniaxial nematic order. 
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FIG. 2. (color online) Binding energy E\, as a function of the 
scattering length as for extreme prolate (EP), extreme oblate 
(EO) and spherical (S) gauge field configurations. In the EP 
case, a bound state is obtained only for l/a^ > 0, while for 
the other two it is obtained for every Os. In the BCS side, 
-Eb depends exponentially on Aa^ for the EO configurations, 
while it is a power law in the spherical case. 



Spherical (S) 



This most symmetric GFC is characterized by A^; — 
Ay ~ ~ ^ ^^"^ marked by the red line (along the 
body diagonal) in Fig.[IJa). Apart from translation and 
time reversal, this GFC has global rotational symmetries 
about all three axes generated by Ji — Li 4- ^r^. Again, 
we find that a two-body bound state appears for any 
scattering length, i.e., a^c = 0^. Also, we obtain a 
closed form expression for the binding energy (referred to 
the scattering threshold Etu — ~\) for any scattering 
length 




(13) 



An interesting aspect of this result is that, for a small 
negative scattering length (BCS side), the bound state 
energy depends algebraically on and A, 



where the subscript r on the second term indicates that 
the spin quantization axis is along r. The wave func- 
tion is made of two pieces. The first piece corresponds 
to a J = state constructed out of L = orbital state 
and a spin singlet, while the second piece is a J = 
state obtained by fusing an L = 1 orbital state and a 
spin triplet state. Furthermore, orbital wave functions 
of both pieces are non-monotonic, i.e., they have spa- 
tial oscillations. This owes to the existence of two length 
scales determined by Ef, and A. While the former dic- 
tates the exponential decay of the wave function, the lat- 
ter determines the period of its spatial oscillation. This 
observation also applies to the wave functions discussed 
above for the extreme prolate and extreme oblate cases. 

It is interesting to note how the bound state evolves as 
we go from the EP to the EO GFC along the path in GFC 
space indicated in Fig.[T|a). The prolate side of the path 
which has a bi-axial nematic order, is separated from the 
oblate side with a uniaxial nematic order by the spherical 
configuration. For the spherical configuration the bound 
state is fully (spatial + spin) rotationally symmetric. 



Eb 



Aa. 
~3~ 



(14) 



D. Generic GFC 



i.e., a deeper bound state than the EP case is obtained 
(see Fig. [2]) in this case. For small positive a^, the leading 
term in the binding energy is that in the free vacuum. 
The bound state is a J singlet and has the wave function 



. Ar 1 I ^ , , 



^ 1\ Ar 



A Ar \ e 

-— cos —- — , 

\/3 y/Ebr 



For a generic GFC, the critical scattering length for 
the formation of a bound state can be expressed as 



XGsc = J^W 



(15) 



where is a dimensionless number- valued function of the 
unit vector A. The function T has to be obtained numer- 
ically. We find that is a non-positive function, i.e., for 
a generic GFC, the critical scattering length is negative, 
i.e., on the BCS side of the resonance. In other words. 



5 



-1 



n 



FIG. 3. (color online) Renormalization group flow diagram 
(schematic) of the two body problem. The point Vp cor- 
respond s to the free fixed point and to resonance in 
vacuum.'^^'^ The non-abelian gauge field is a relevant op- 
erator at these fixed points as indicated. 



the strength of the critical attraction required to produce 
a two body bound state is reduced by the presence of a 
generic gauge field. Fig. [l|b) shows the evolution of asc 
along the great circle connecting the EP state to EO state 
for a fixed gauge coupling, illustrating that prolate GFC 
has a negative asc, while any oblate GFC has a vanish- 
ing Use- In summary, the two-body boundstate appears 
at resonance (l/a^c = 0) for EP GFCs marked by the 
blue lines (along the axes) in Fig.jlja). For spherical (S) 
GFC marked by the red line (along the body diagonal) 
and for oblate GFCs marked by the planes bounded by 
the green (including EO) lines (along 45° to the axes on 
the coordinate planes) and the red line (along the body 
diagonal), a^c vanishes, i.e., any attractive interaction, 
however small, will force a bound state for the two body 
problem. 



IV. QUALITATIVE DISCUSSION 

We now discuss the physics behind these results. In 
the free vacuum, a renormalization group analysis of the 
field theor y of the two body problem with the contact 
interactiorf22l2i] provides two fixed points. The first is a 
stable one at f ^ = describing two free fermions, and 
the second, an unstable one u^j = — 1 (in suitably cho- 
sen units) corresponding to the resonance; see fig. [s] In 
free vacuum, a contact interaction parameter v near Vp 
fiows toward Vp and hence has similar physics as two 
free fermions. This corresponds to the fact that suffi- 
ciently strong attraction is required {v < —1) to produce 
a bound state. Consider now a situation with a v near 
Vp and a non-vanishing spherical gauge field with a cou- 
pling strength A. We see immediately that Rashba term 
described by the coupling A is a relevant operator and 
the flow takes the system away from the free fixed point 
(see Fig. [s]) suggesting that even a small A has a drastic 
effect on a system near the free fixed point. 

A deeper understanding can be obtained by consider- 
ing the density of states g{e) of a single fermion mov- 
ing in a gauge field, since it determines the density of 



states of two non-interacting fermions with zero center 
of mass momentum. One can be easily obtain analytical 
expressions for the density of states for the high symme- 
try GFCs. The gist of those formulae is that near the 
scattering threshold, for the high symmetry GFCs dis- 
cussed above 



< A(constant) 

Its 



for EP 
for EO 
for S 



(16) 



In all three cases g{e) — ^/e as e — oo. It is therefore 
clear that the infrared behavior of the density of states 
is behind the results presented hitherto. This motivates 
us to construct a model with density of states given by 



9{e) = 



/2e 



e < eo, 
if e > eo. 



(17) 



where O is the unit step function, 7 is an exponent that 
determines the infrared behavior of g, and Eq is an en- 
ergy scale (crudely equal to A^) at which the density of 
states is restored to that in the free vacuum. Note that 
7 — |,0, — ^ qualitatively reproduces, respectively, the 
density of states corresponding to EP, EO and S GFCs. 
We can readily calculate the critical scattering length 
that obtains a bound state as 



27-1 



9(7) 



(18) 



We see immediately that a^c vanishes whenever 7 is non- 
positive as is the case for the EO and S configurations. 
For the EP configuration, the critical scattering length 
asc —00 consistent with the results presented earlier. 
For a generic GFC, the infrared density of states has 
a narrow y^-regime, followed by a regime with nearly 
constant density of states - this can be modeled in this 
simple picture using a 7 that satisfies < 7 < ^, the 
precise value of 7 being dependent on the direction A 
in the GFC space. We find that Usc is negative, again, 
consistent with our calculations. 

This simple analysis allows us to uncover the physics 
behind the phase diagram of Fig. [l] Highly symmet- 
ric GFCs drastically modify the low energy density of 
states owing to the degeneracies induced in the resulting 
one particle levels. It is in this sense that the Rashba 
term is a relevant operator as mentioned in the discus- 
sion above. For highly symmetric GFCs, the enhanced 
density of states at low energies strongly promotes bound 
state formation in the presence of an attractive interac- 
tion. 

The particular type of nematic spin symmetry in the 
bound state arises so as to optimize the kinetic energy. 
The spin-orbit interaction mixes the singlet and triplet 
sectors of the two particle system, and the particular ne- 
matic symmetry obtained in the bound state enables the 
orbital wave function to sufficiently "sample" the attrac- 
tive interaction at a minimal cost in kinetic energy. 
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V. EXPERIMENTAL DIRECTIONS AND 
OUTLOOK 

Our predictions can be readily tested by experiment. 
A clear signature of the bound state formation can be 
obtained from the measurement of en erg^i^ of the gas at 



high temperatures. As suggested by Ho and MuelleijP , 
a large value of the second virial coefficient which char- 
acterizes the interaction energy, is obtained for an inter- 
acting Fermi gas in free vacuum near resonance. Our 
results suggest that in presence of a generic gauge field, 
such large corrections to energy can be observed on the 
BCS side i.e., for negative scattering lengths near a^c at 
the onset of the bound state. The quantitative value of 
age will be affected also by the attraction in the triplet 
channel, but our predictions can be tested qualitatively. 

Our results suggest that the BCS-BEC crossover in the 
presence of a non-Abelian gauge field will be drastically 
altered and in particular the "crossover regime" should 
shift to the BCS side for a generic GFC. There are many 
other novel effects of the non-Abelian gauge field in the 
many body context such as transition in the topology of 
the Fermi surface with increasing fiUin^SU^ Moreover, the 



two-body bound-state wave functions provide a clue to 
the nature of the Cooper pair wave functions. Clearly, 
this will lead to superfluild states with interesting pair- 
ing wave functions (such as the extended BW and ABM 
states found here with associated nematic orders) and 
concomitant excitations. 

The simple model we have presented in Sec. |IV| sug- 
gests that our conclusions will also apply to systems with 
a larger gauge group (such as SU(7V)). Investigations 
along these lines should lead to interesting new possibil- 
ities with cold atom systems. 
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